Mean field solutions of kinetic exchange opinion models 
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We present here the exact solution of an infinite range, discrete, opinion formation model. The 
model shows an active-absorbing phase transition, similar to that numerically found in its recently 
proposed continuous version (Lallouache et al, Phys. Rev E 82, 056112 (2010)). Apart from the 
two-agent interactions here we also report the effect of having three agent interactions. The phase 
diagram has a continuous transition line (two agent interaction dominated) and a discontinuous 
transition line (three agent interaction dominated) separated by a tricritical point. 
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I. INTRODUCTION 

The dynamics of opinion formation in a society and emer- 
gence of consensus arising out of cooperative interactions 
between the agents are being extensively studied recently 
Due to its cooperative nature, statistical mechan- 
ics provides the required tools to study such systems. 
Although many intricacies of real societies are lost in the 
process, such minimal modeling often yields essential fea- 
tures in terms of their social as well as physical aspects. 

The key feature in modeling opinion formation is to 
quantify opinions in terms of real numbers. Depending 
upon context, opinion is often quantified as discrete or 
continuous variables between two or more choices. Also 
the process of interaction between the agents is a vital 
ingredient. While several choices to model such inter- 
actions exist, one way is to consider an interaction as a 
'scattering process' where the agents are stochastically 
influenced by each other's opinions (see, for example, @- 

m 

Recently an opinion formation model [13j based on 
such a 'kinetic exchange process' between two individuals 
was proposed [the Lallouache-Chakrabarti-Chakraborti- 
Chakrabarti (LCCC) model hereafter]. Resembling the 
model for wealth exchange in a society 14j, this model 
has a single parameter that determines the 'conviction' 
of an individual. It was shown that beyond a certain 
value of this conviction parameter, the society reaches 
a 'consensus', where one of the two choices (positive or 
negative) provided to the individuals prevails, thereby 
spontaneously breaking a discrete symmetry. 

A generalization was proposed subsequently [l5j . in 
which the 'self-conviction' and the ability to influence 
others were taken as independent variables. This two- 
parameter model has a simple phase boundary along 
which apparent non-universality was observed (for de- 
tailed discussion on the critical behavior of a class of 
model of this kind see [TH, H3])- Subsequent extensions 
in terms of including "noise" [l8| and also studying a 
generalized map version [l9[ for this class of models were 
done. 

In the present paper, a discretised version of the LCCC 
model is exactly solved. As this is an infinite range 



model, the solution is essentially mean field. Also a 
generalisation in terms of three-agent interactions is pro- 
posed. From the expression of the order parameter it is 
seen that for pure three-agent interactions (Model I) the 
transition is discontinuous (giving hysteresis like behav- 
ior as well) but for mixture of two-agent and three-agent 
interactions (Model II) , the nature of transition depends 
on the relative probabilities of the two types (two-agents 
and three-agents) of interactions. 

As we will see, the present version of the model is 
a (three-state) probabilistic cellular automaton (PCA) 
(see, for example, (2fj| - [24l ]). Such types of models have 
been used in sociophysics for a long time starting from 
Schelling [25| (see also [26l - l29j to name a few) . 

The paper is organized as follows: In Sec. II the exist- 
ing (LCCC) and the proposed models will be described. 
In Sec. Ill the mean-field calculations to find the ex- 
pression for order parameter is presented. Then in Sec. 
IV the three-agent generalisation is introduced with the 
analysis of the order of transition. In Sec. V, we discuss 
the phase boundary obtained for the model with both 
two-agent and three-agent interactions. Finally, the re- 
sults are discussed. 



II. KINETIC EXCHANGE OPINION MODEL 
AND ITS DISCRETE VERSION 

Consider a simple pair-wise interaction model for opinion 
formation (LCCC) in a well-connected group of individ- 
uals. The opinion of an individual in the society is rep- 
resented by a real number which can continuously vary 
within the limit — 1 < 0$ < +1. At any time t an agent 
with opinion Oi(t) interacts with another randomly cho- 
sen agent with opinion Oj(t). After the interaction the 
i-th agent retains a fraction of her own opinion (which 
depends on the agent's 'conviction') and is stochastically 
influenced by the jth agent. The amount of the influence, 
of course, depends upon the jth agent's 'conviction'. The 
dynamics of the LCCC model evolves following the equa- 
tion 



o.it + 1) = AjOj(t) + Xj€Oj(t). 



(1) 
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The parameter A m represents the conviction of m-th 
agent and e is a stochastic variable uniformly distributed 
between [0,1]. If the opinion of an agent reached the 
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higher (lower) extreme +1 ( — 1), then of course its opin- 
ion value was prevented from further increase (decrease). 
N such exchanges (where N is the total number of 
agents) constitute a single Monte Carlo (MC) time step. 
For simplicity, the society was assumed to be homoge- 
neous in the sense that all A m 's were same (say, A). 
(Note that in Ref.[l3j], simultaneous changes of opinions 
for both the ith and jth agents were considered, whereas 
Eq. ([T]) describes a single change.) As long as the dy- 
namics is of infinite range (any agent can interact with 
any other), this difference does not affect the results, ex- 
cept for the fact that the time is rescaled (2N updates 
for a MC step). 

The steady-state characterization of this model was 
done using two measures. One is the average opinion 

JV 

of the agents O = 0j representing the measure of 

i=i 

global consensus and the other is the fraction of agents 
having extreme opinions 

C = /i + /-i, (2) 

where /i and /_i are the fraction of agents having opin- 
ion value +1 and —1 respectively. 

Extensive Monte Carlo (MC) study yields that in 
the LCCC model, above A c ~ |, O ^ and below A c , 
O = 0. Furthermore, in the disordered phase, which is 
rather special here, o; = for all i. This can be con- 
sidered as a neutral state of the agents (avoiding such 
neutral states are possible in some cases [IH). Similar 
behavior was also obtained for C. As in usual critical 
phenomena, the relaxation time shows divergence from 
both sides of the critical point following a power law, 
having same exponent value on either side of the criti- 
cality. Although nothing could be predicted about the 
critical exponent values, a mean-field like analysis gave 
A c = 2/3 for the LCCC model (for detailed discussions 
see Ref.[I3). 

In its generalization [l5j], it was argued that the 'self- 
conviction' A of an agent need not, in general, be equal to 
her ability to 'influence' others. In its generalized form, 
therefore, the dynamical exchange process reads 

0i (t+ 1) = Ao;(t) +fJ£Oj(t), (3) 

where /i represents the j-th agents ability to influence 
others. In the limit A = /i one recovers the LCCC model. 
For this generalised model, there is a phase boundary 
in A — [i plane, having the equation A c = 1 — /i c /2. The 
values of the critical exponents along this phase boundary 
were reported to have strongly non-universal behavior 
[l5j for O and weakly non-universal behavior for C. 

The above-mentioned models defy simple analytical 
treatments to find the order parameter as long as Oj's 
are continuous. But it is often the case in a society that 
the opinions can take only discrete values (e.g., voting 
'yes' or 'no' for a referendum, or voting in a two-party 
political scenario etc.). While retaining the social inter- 
pretation, it significantly simplifies analytical treatment. 



To that end the following modifications are made in the 
present study: The dynamical exchange equation [Eq.([T])] 
remains the same, but we make A stochastic in the sense 
that we put A = 1 with probability p and with proba- 
bility 1 — p. Also, the parameter e is either 1 or with 
equal probability. This is a discretized version of the 
LCCC model where the agents have three possible opin- 
ion values (i.e., Oi £ { — 1, 0, +l}Vi). This is, therefore, a 
three-state probabilistic cellular automaton (PCA). 

The above mentioned modifications lift the 'homoge- 
neous society' (all agents having same 'conviction') as- 
sumption and on the other hand keep the opinion val- 
ues discrctiscd. However, the inhoinogeneity is the sim- 
plest of its kind: only 'high' and 'low' convictions are 
present. The agents can change between these two states 
randomly in time (A is annealed variable). The case of 
quenched A, in this case, is a trivial limit where order 
parameter becomes simply proportional to p. 

In the case of the generalised version [Eq.Q] the ad- 
ditional change is that like A, we put fi — 1 with prob- 
ability q and with probability 1 — q. However, regard- 
ing its variation in time, it is explicitly checked numeri- 
cally throughout the paper that the results do not change 
whether /i depends on time or not [at least in the mean- 
field (MF) limit] . Therefore, to facilitate analytical treat- 
ment, it is assumed to be randomly varying with time. 

Finally note that the relevant parameters of the prob- 
lem here will be p and/or q, which essentially specify the 
average values of A and fi respectively. Also, the present 
versions are three-state PCA, but one can prove (see Ap- 
pendix A) that in the steady state opinions of opposite 
signs do not coexist (both in 'space' and time). There- 
fore, all the subsequent results will hold for a two-state 
PCA as well. 



III. MEAN FIELD SOLUTION OF THE 
DISCRETE MODEL 

It is our intention to find an expression for the order 
parameter O in terms of p and to find out the order 
parameter exponent (3 defined as O ~ (p — p c )^ . Subse- 
quently a similar expression in terms of both p and q for 
the generalised model (see Eq.Q) is also presented. 

Let fo, fi and /_i be the fractions of agents having 
opinions 0, +1 and —1 respectively. Now, since the inter- 
actions are only pair- wise and both A and fi can take only 
two values, one can enumerate all possible interactions 
between all possible pairs, which contribute to increase 
and decrease of the order parameter. For example, the 
probability that both the agents in an exchange process 
have opinion +1 is ff. Now the probability with which 
one agent shifts her opinion to is (1 —p). Therefore the 
process (1,1) — » (0,1) has probability ff(l — p). This 
process, of course, contributes to decreasing of the order 
parameter. One can enumerate all the eight processes 
that contribute to changing the average order parameter 
(O = fx — f-i) and can therefore write the evolution 
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FIG. 1. (Color online) Variation of the order parameter O 
for discrete LCCC model and its generalised version (in the 
limit p — q). Analytical expressions given by Eqs. ([8| (lower 
curve) and (j 10|) (upper curve) are in good agreement with the 
simulation points. The critical or symmetry breaking point 
(SBP) p c = 2/3 is indicated. For the simulations, N = 10 5 
agents were considered, The initial conditions are ordered (all 
+ 1 or all —1). The results are time averaged (100 MC steps) 
after the system reaches a steady states (roughly after 10 3 
initial MC steps) and about 100 ensemble averages (different 
starting configurations) were taken. 
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FIG. 2. (Color online) Variation of the order parameter O 
for pure three agent interactions. The solid line is given by 
Eq. (fT6)> i.e., the r = 1 limit of Eq. I|22[). In the simulation 
points, when the initial condition is an ordered state (all +1 
or all —1), the discontinuous jump occurs at p c i — 8/9, while 
starting from the random initial condition, the jump occurs at 
p C 2 ~ 0.930 ± 0.005. This clearly shows a hysteresis behavior 
(or, initial condition dependent behavior), which signifies the 
'resistance' of the society to global changes in the opinion. 
N = 10 5 agents are considered for the simulation points. 



equation as 



f = £x(l - p) + f-iMl - | ) + + /-i/od - P) 

~/i 2 (l - P) - A/-i(l - |) - /i/o(l - p) (4) 

In the steady state, the left hand side will be zero. Then 
after simplification one arrives at 

A 2 (l-p)+/o/-i|+/o/i(l-p) 

= /_ 1 (1-P) + /o/i| + /o/-i(1-p) (5) 

This gives either fx = f—i, (which implies disorder) or 



/o 



2(1 -p) 



(6) 



It is possible in this case to show explicitly (see Ap- 
pendix: A) that in the ordered state = (making 
the steady state value of O and C identical in this and 
for all subsequent discussions also). This condition and 
the disordered state condition (fx = should both 
be valid at the critical point. This is possible only when 
fx = /_i = at the critical point. This implies, at the 
critical point fa = 1. Furthermore, for the sake of con- 
tinuity of fx and f—i, fa = 1 for the entire disordered 
phase. This condition along with Eq.([5]) gives p c = 2/3. 

Therefore, the order parameter should be (using fx + 
/-i+/o = l) 



O = ±(1 - fa) 



(7) 



where the sign will depend on whether fx or /_i is non- 
zero in the ordered (symmetry-broken) phase. Using 
Eq. the above expression yields 



o = ± 



3(p 



(8) 



Therefore, Eq.© gives /3 = 1 (since p c = 2/3). 

This can of course be generalized for Eq.©. Straight- 
forward algebra would yield 



/o 



(p-l)(?-2) 



pq 



(9) 



As before, in the disordered phase /o = 1, which yields 
the equation for the phase boundary in the p-q plane as 
p c = 1— • This gives the expression for order parameter 
as 



= ± 



2(p-p c ) + - Qc) 



pq 



(10) 



Therefore, no matter through which path and which 
point the phase boundary is crossed, the order param- 
eter exponent is /3 = 1. The discretized version of the 
LCCC model presented here has a transition, which is 
mean-field active-absorbing type [23|, [24j]. Accordingly 
/3 = 1 is obtained. 

Of course we do not expect to get Eq.® from Eq.(fT0|) 
by putting p = q, as this would only mean (A) = (/i) and 
not A = /i. 

The above-mentioned results arc verified using numer- 
ical simulations. Ordered initial conditions (all +1 or all 
— 1) were taken. The system was allowed to relax un- 
til steady state was reached (checked by monitoring the 
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time variation of the order parameter). For N — 10 5 , 10 3 
initial time steps were discarded for relaxation (depend- 
ing on the proximity of the model parameters p, q with 
their critical values, this has to be increased upto 10 4 for 
some of the points presented). Then a time average was 
taken for about 100 MC steps (1 MC step being N update 
attempts). The whole process was averaged over 100 en- 
sembles. Note that a random initial condition (+1, —1, 
with equal density) will give same steady-state values for 
the order parameter. But the relaxation time will be 
higher. 

In Fig. [TJ Eq. © is compared with Monte Carlo sim- 
ulations to find good agreement. 



IV. MODEL I: BEYOND PAIR- WISE 
INTERACTIONS: THREE-AGENT 
INTERACTION AND DISCONTINUOUS 
TRANSITION 

In all previous studies [l3|, UH, [TtJ regarding the kinetic 
exchange processes mentioned here, interactions were al- 
ways taken to be occurring between two agents. This 
is partly because two-body exchange is the simplest and 
also because in the energy exchange of ideal gas only two- 
body interactions are important. But in opinion forma- 
tion, exchange between more than two agents is perfectly 
possible. So we intend to investigate the effect of such 
interactions in opinion formation. 

The simplest possible generalisation toward many- 
body interaction is to consider three-body exchange. In 
doing so, the following strategy is followed. Three agents 
are chosen randomly. Then one agent modifies her opin- 
ion according to that of the other two only when the 
other two agrees among themselves. If they do not agree 
the first agent considers the group to be 'neutral'. Math- 
ematically this can be represented as 



0i (t+l) = Xoi(t) + Xe9 jk (t), 



where, 



(*) 




if Oj(t) = o k (t) 
otherwise. 



(11) 



(12) 



Before proceeding further, it is to be noted that 9jk if) 
takes the value +1,-1 and with probabilities ff, f^_ 1 
and 1 — (ff + /^j) respectively. Then just as Eq.([T|) was 
treated, one can enumerate all exchange processes that 
contribute to increase and decrease in the order parame- 
ter. The time evolution equation is 

■ ¥ =/- 1 (1-p) + /-i/i(1-2) + — 

+/_! [i - (A 2 + A)] (i -p) - ff(i -p) - /i/ 2 i(i - |) 



fof-iP 



/l[l-(/l 2 + £l)] (1-P) 




200 



FIG. 3. (Color online) The time dependence of O and C 
are plotted for two different values of p (0.90, 0.92) both in 
the range p c i > p > p C 2 for random (o, = 1 and —1 with 
equal probability) and fully ordered (oi = 1 for all i) initial 
conditions in the pure three agent interaction case. O and 
C are identical by definition for the ordered initial condition. 
For the random initial condition, it is seen that both C and O 
go to zero exponentially (linear in log-linear scale) confirming 
the initial condition dependent behavior in this region of p 
values. Number of agents is N = 10 s in the simulations. 



Again, in the steady state left hand side should be zero, 



giving 



/i 3 (i-p) + /i/Mi 



f f 2 p 

JoJ-i 2 



+/i[l-(A 2 + £i)] (1-P) 
= / 3 1 (l- P ) + /_ l/l 2 (l-| 

f-i [i-(A 2 + /-i)] (1-p). 



/o/i 



IP 



(14) 



This gives cither /i = /-i (which implies disorder) or in 
the ordered state 



= 1 _3yV£8/9 
/0 2 2^p ' 



(15) 



where we have used fif—i = as a simplifying assump- 
tion which is numerically verified here. We have ne- 
glected one solution of /o in which it increases with re- 
spect to p in the ordered phase (because in the limit 
p = 1, we expect /□ = 0). Using this, the order parame- 
ter takes the form 



= ± 



3^p^8/9 

2Vp , 



(16) 



(13) 



Clearly, the above equation gives real values for O only 
when p > 8/9. Therefore, for p < 8/9 the only real 
solution can be /i = /_i i.e., O = 0. But from the 
form of Eq. ([T(3f it is clear that in the ordered phase, the 
minimum value of O can be 1/2. Therefore, the order- 
disorder transition is necessarily discontinuous. 
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FIG. 4. (Color online) The phase diagram of the model with 
mixture of three agent and two agent interactions in the p-r 
plane, where r denotes the fraction of three agent interac- 
tions. The continuous transition line follows Eq. (|23|) . while 
the discontinuous transition line follows Eq. (|24[l . Clearly, in 
the limits r = and 1 the transition points are 2/3 and 8/9 
respectively, as is expected from the discussions in the text. 



To verify this and to see the initial condition dependent 
behavior, MC simulations were performed. Depending on 
the initial condition, the discontinuous jump from order 
to disorder happen at two different points, thus showing 
hysteresis behavior (see Figj2j . When the initial condi- 
tion is ordered, Eq. (fl"6|) is followed upto p c \ = 8/9. After 
that the order parameter jumps to zero. On the other 
hand, when the initial condition is random (having al- 
most equal number of agents having opinions of opposite 
signs) then O = upto p c2 ~ 0.930 ±0.005 and then sud- 
denly jumps to the ordered (symmetry broken) phase. Of 
course p C 2 is the symmetry breaking point. Note that the 
estimation oip c i is entirely numerical here. Fig. [3] shows 
the time dependence of O and C in the region where 
Pd > p > Pci. It is seen that for random initial condi- 
tion, both C and O go to zero. This signifies that Oj = 
for all % in this case, while those curves for the ordered 
initial condition shows saturation to a non-zero value. 
This clearly confirms the initial condition dependent be- 
havior in this region of p values. For one particular point 
(p = 0.91) this is checked for large enough size (N = 10 6 ) 
for two different initial conditions. It is clearly seen then 
the long time saturation values are very different, as is 
expected. 

This initial-condition-dependent behavior can be 
somewhat understood as follows: When the initial con- 
dition is random, the two agents chosen for the inter- 
action with the first one have rather smaller chance of 
agreement among themselves. So the first agent sees 
the group as neutral (i.e., with opinion value zero). In 
this way, the agent finds itself in an environment that 
effectively has a lot of neutral agents. Therefore, un- 
less the agent's conviction is very high (turns out to be 



FIG. 5. (Color online) Variations of the order parameter for 
r = 0.2(< rt), r = 0.5(= n) and r = 0.8(> n). The contin- 
uous lines are analytical results (see Eq.([22}) and the points 
are simulation results with N = 10 5 . For the first two curves, 
continuous transitions are seen, as expected. For the last one 
(r = 0.8) a discontinuous transition with signature of hystere- 
sis is seen. Inset shows the log-log plots of O versus p—p c near 
the transition points for the continuous transitions. From the 
slopes of the curves the order parameter exponent /3 is found 
to be 1 and 1/2 for r = 0.2 and r = 0.5 respectively. 



about 0.93), all of them succumb to this neutral environ- 
ment and an absorbing phase is reached. The situation, 
however, is very different when the initial condition is 
ordered. Here the environment is rather supportive in 
maintaining the prevailing opinion in the society (as is 
generally seen). Therefore, intuitively we expect that 
only at a lower value of conviction, the prevailing opin- 
ion will decay. This gives rise to the initial-condition- 
dependent behavior. 

Initial-condition-dependent behavior in opinion forma- 
tion has been studied before in different contexts. For 
example, in Ref. [29j, 'hysteresis' was observed while 
modeling the influence of a strong leader in the society. 
In general, these behaviors signify the resistance offered 
by the society to changes in the global opinion (both from 
order as well as disordered state), despite the fact that 
the very reason for its formation has lost its relevance. In 
the present case too, the 'hysteresis' (this word here does 
not have the usual meaning, say, that in the magnetic 
system) loop area is somewhat a measure of this social 
tolerance and resistance. 

A similar approach can be made for the generalised 
case described by Eq.©. However, instead of doing so, 
one could also look at another limit of Eq. © , where q = 1 
(studied in Ref. [l7[ as Model C). One can show that even 
in this limit, a discontinuous transition can be obtained 
with p c i — (2 + v2)/4. So one would in general expect a 
discontinuous transition for all ranges of Eq. © . 
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V. MODEL II: MIXTURE OF TWO-AGENT 
AND THREE-AGENT INTERACTIONS: PHASE 
DIAGRAM AND TRI-CRITICAL POINT 

Let us now discuss how robust is this discontinuous tran- 
sition. In the above analysis it was assumed that only 
three-agent interactions are present as opposed to the 
previous cases, where only two-agent interactions were 
considered. Here we consider a situation where both 
two-agent and three-agent interactions are allowed. In 
principle interactions of all sizes should be allowed, but 
this is the simplest generalisation one could make. 

With probability r an exchange process is three-agent 
and otherwise it is two- agent. The exchange equation is 
the same as Eq. fTTj) but now clearly 

!1 with probability rfl + (1 —r)fi 
-1 with probability rf\ + (1 - r)/_i 
otherwise. 

(17) 

With this, one may enumerate all possibilities of increase 
and decrease of the order parameter and the time evolu- 
tion equation will be 

^=/-i [rf\ + (1 - r)f-x] (1-p) 

+h [rf% + (1 - r)/_x] (1 - | ) + f -f [rfl + (1 - r) h ] 

[1 - (fl + f- x )r - (1 - r){h + f-i)} (1 - p) 
-h Vfi + (1 - r)h] (1-P) 

/— i [rfl + (1 - r)/i] (1 - f ) - f ~f [rfl, + (1 - r)/_x] 

-fi [1 - (fl + f\)r (1 - r){f x + f-i)] (1 - P). (18) 

In the steady state, the left hand side will be zero, giving 
(assuming f\f-\ — 0, which can be verified numerically) 

h [rfl + (1 - r)h] (l-p) + f -f [rf\ + (1 - r)/_i] 
+/i [1 - (fl + f\)r - (1 - r)(/i + (1 - p) 

= S-x Vf\ + (1 - r)/_ x ] (1 -p) + ^ [r/f + (1 - r)/ x ] 
[1 - (fl + f-i)r - (1 - r)(/x + /-0] (1 - p). (19) 

On simplification, this yields either /i = /— i which im- 
plies disorder, or 



pr, 2 P/o , , n 

y/o- — + i-p = o, 

which gives (the only relevant solution) 



_ 1 y/p 2 /4-2pr(l-p) 

JO — 7, ' 

Zr pr 



(20) 



(21) 



Again as before 



= ±| 2r-l + ^74EWEg . (22) 



The sign is determined from the fact whether /i or /_ i is 
non-zero in the ordered phase. However, the terms inside 
the bracket essentially represent a probability (1 — fo) and 
therefore, should be a positive and real number. Now, 
when r < 1/2 for the quantity to become positive, p 
has to be sufficiently high (> p c , say). Before that the 
disordered solution (O — 0) remains valid. Therefore, 
the ordered-disordered phase is separated by a line in 
the p — r plane, which is obtained by setting the quantity 
inside the bracket to zero. That critical line is 



Pc 



for 



1 



(23) 



where p c , r c are the set of values satisfying the above 
equation where a continuous transition takes place. Note 
that the above condition also keep O real, which is the 
other requirement, as long as r < 1/2. In Eq.(f22j) one 
can put p = 2/(3 — r) + A (where A — > 0) and show 
that the leading order term comes out to be linear in 
A, implying /3 = 1 along this line. This critical line, of 



course, terminates at (p t 



When r > 1/2, the ordered-state solution Eq. (j2"2")l 
can be valid whenever it gives real values for O. For this, 
the phase boundary is obtained by setting the quantity 
inside the square-root to zero. The said phase boundary 
is 



Pc 



8r c 



l + 8r c 



for r > 



1 



(24) 



When r > 1/2, the minimum value possible for O from 
the ordered state solution (Eq. (|2"2"|) ) is greater then zero. 
Therefore transition across this line is necessarily discon- 
tinuous. This discontinuous nature is verified numeri- 
cally. The amount of discontinuity, of course, is given 
by 1 — gr, which is maximum (1/2) for pure three agent 
interactions (r = 1). A 'hysteresis' like behavior, as dis- 
cussed in the preceding section, is also seen. Note that 
the phase boundary equations correctly give p c = 2/3 
and p c — 8/9 limits respectively for r = (from Eq.([23"|)) 
and r = 1 (from Eq. pjjl ). 



The point (p t 



i) is special where the critical 



line terminates. It is a Tricritical Point (TCP). As is seen 
generally, at the TCP the exponent values are different. 
Clearly, 



O ~ VP - 4/5 



(25) 



giving fir cp — 1/2, which is different from fi = 1 found 
along the critical line. 

The phase boundaries and the tri-critical point are 
shown in Fig. |4j In Fig. O the order parameters are 
plotted for three different values of r, one on either side 
of the TCP and one at TCP. The inset shows the differ- 
ent values for the exponent j3 when the r < r t (/3 = 1) 
and when r = rt ({3 — 1/2). 

To find the other exponent values that characterize this 
TCP, one can use the following scaling relation for the 
order parameter 



2r 



pr 



0(t) « t- s F (V^iiA,^ 2 /^) 



(26) 
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FIG. 6. (Color online) Data collapse for finding i/u (see 
Eq.p6]l) for different p values for r = 0.5 (p c = 0.5). The 
estimate is v\\ — 1.00 ± 0.01. Inset shows the uncollapsed 
data. TV = 10 s agents are considered, initial condition is or- 
dered (all +1 or all -1). 
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FIG. 7. (Color online) Data collapse for finding z for different 
system sizes (TV = 500, 1000, 5000, 10000, 50000) at p = 0.5 
and r = 0.5 (TCP). The estimate is z/d = 0.666 ± 0.001. 
TV = 10 5 agents are considered, initial condition is ordered 
(all +1 or all —1). Inset shows the uncollapsed data. 



where A = p—p c , v\\ is the time-correlation exponent, z is 
the dynamical exponent and d is the space dimension. At 
the critical point, the order parameter follows a power- 
law relaxation 0(t) ~ t~ s (see inset of Fig. [5] ) with 
S = 0.50 ±0.01. 

One could then plot 0{t)t against t(p — p c Y 11 ■ By 
knowing <5, vu can be tuned to find data collapse. From 
Fig. [5] the estimate of v\\ is 1.00 ± 0.01. Similarly, one 
can plot 0(t)t s against t/N z ^ d . Again by tuning z, data 
collapse is found (see Fig. [7]). The estimate of z/d is 
0.666 ± 0.001. To find z one should put d = 4, which is 
the upper critical dimension. This gives z ~ 8/3. Similar 



analysis for r < i gives z/d 1/2, here also by putting 
d = 4 one gets z « 2, which is expected for directed 
percolation (DP). 

The scaling relation 6 — fi/v\ is clearly satisfied here. 

Therefore we see that the nature of the transition is ac- 
tually determined by the relative probabilities of the two- 
agent and three-agent interactions. Also the exponent 
values at the tri-critical point are different from those 
along the critical line. 



VI. DISCUSSION 

In the first part of this paper an exact solution of the 
discretised version of a recently proposed model (LCCC) 
[lH for opinion dynamics is given. The model being infi- 
nite ranged, the solution is essentially mean field. The in- 
teractions, as in its simplest form, are between two agents 
[see Eq.([T])]. The exchange process is such that an agent 
has a 'conviction' with which she retains her opinion and 
also gets influenced (stochastically, because it is other- 
wise impossible to incorporate all social complexities in- 
volved in such processes) by the opinion of one randomly 
chosen agent. It was shown [l3|, [l5[ from extensive MC 
study that beyond a certain value of the 'conviction' pa- 
rameter the society undergoes a phase transition from 
disordered to ordered state (where consensus is formed) . 
In the present study that behavior is shown analytically 
[see Eq.©] for a discretized version of the model in mean- 
field limit (which is exact here). The order-parameter 
exponent has been found to be (3 — 1. Even for the gen- 
eralised version [l5[ [see Eq.Q] this exponent remains 
the same along the phase boundary (belonging to DP 
universality class). 

Thereafter, a generalisation of this model for three- 
agent interaction is reported (Model I). There is, of 
course, no single choice for this kind of generalisation. 
But here we have taken a plausible strategy in which an 
agent can be influenced by the opinions of two other ran- 
domly chosen agents only when those two agents agree 
among themselves (have the same opinion) otherwise the 
first agent considers the group as 'neutral' [see Eq. fTTj) ]. 
This generalisation has led to an interesting behavior in 
terms of the order of the transition. It is seen if all in- 
teractions are three- agent, a discontinuous transition is 
obtained [see Eq. (|16[l ] and an initial condition dependent 
behavior was also observed (Figj2]). It is to be noted that 
'hysteresis'-like behavior in opinion models have been re- 
ported before in other contexts (see e.g., Refs. (29l. l3fj|). 
In general, this behavior in opinion dynamics models can 
be taken as a signature of the tolerance of the society, or 
in other words, its resistance to changes in global opinion 
(as is also indicated in Ref. 29] ) . Although a direct cor- 
respondence to a measurable quantity cannot be made 
from these simplified models, qualitatively this loop-area 
is somewhat a 'measure' of this social 'tolerance' men- 
tioned above. 

It is important to find out how much of this discontin- 
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uous nature is generic or it is an artifact of the restriction 
of only three-agent interaction, as invoked by Eq. (jTTJ) . 
Of course it is not possible (or at least it is very diffi- 
cult) to allow interactions of all sizes as it should be in 
a real society. But at the very least one can allow both 
two-agent and three-agent interactions with some proba- 
bilities (Model II). In doing so it is found that up to the 
point when the probability of three-agent interactions is 
below 1/2, the transition is continuous [phase boundary 
given by Eq. (|2"3"jl ] and beyond that the transition is dis- 
continuous, [phase boundary is given by Eq.([53]. The 
point where the two-agent and three-agent interactions 
are equally probable, is a special point, because it is 
a tri-critical point. The transition here is continuous, 
but the values of the exponents are different from those 
along the critical line. Along the critical line, the ex- 
ponents are of course mean-field active-absorbing type 
(p = 1, z = 2, v\\ = 1, 6 = 1) but at the TCP they 
are f3 = 1/2, z/d m .666 ± 0.001 2/3 (z w 8/3), 
5 = 0.50 ± 0.01 « 1/2, v\\ fa 1.00 ± 0.01. 

At this point it is appropriate to mention that the 'dis- 
ordered' phase in all these versions is quite special in the 
sense that all the agents are neutral in this phase. One 
can avoid such situation (see Ref. [lH) and make the 'dis- 
ordered phase' have coexistence of opinions of different 
signs. But here such generalisations were not discussed. 
Even without such complexities, which one can add any- 
way with this, one finds intriguing features in this model. 

One may note that while attempting to interpolate be- 
tween a continuous and discontinuous transition, a tri- 
critical point was obtained also in Ref. [3l|. There too, 
the tricritical point was situated at the point where the 
phase boundary changed its curvature. 

Finally, one must also note that in all the above cases, 
the existence of a phase transition has not been proved. 
The two solutions representing ordered and disordered 
states were analytically obtained. The existence of the 
phase transition, for all practical purposes, has been as- 
sumed here. 

To conclude, a mean-field solution for a kinetic ex- 
change model of opinion formation and its phase tran- 
sition in terms of forming global consensus is presented 
here. A genrelization to also include the three-agent in- 
teractions is proposed. Surprisingly, the nature of the 
transition depends on the relative probabilities of the 
two-agent and three- agent interactions. 



Appendix A: Derivation of the polarization 
condition fif—i =0 

Here we show that for the models with two-agent inter- 
actions studied here, in the steady state of the ordered 
phase agents with opinion of opposite signs do not coex- 
ist in both 'space' (all agents) and time. In other words, 
we show here that Jfj v _ x = for any x,y ^ 0. 

For the model defined by Eq. (pTJ) , the time evolution 



equations governing f\ and /_i are 

df-i fof-iP , 2 n \ r f ft P\ 
— = ^--/_ 1 (l-p)-/_ 1 / 1 (l-~) 



dfi fofiP f 2 M \ f * n P\ 



(Al) 



-A/o(l-p). (A2) 



In the steady state we should have 

fof-lP ,2/-, \ f r f-, P\ 

— 2 /-i(l-p) -/_i/i(l--) 



/-i/o(l-p), (A3) 



^~fl(i-p)- A/-x(i-|) 

= /i/o(l-p). (A4) 

Adding the above two equations and imposing the con- 
dition for order state (/i =^ we get 

(1 - p)(l - fo f - 2/i/_x(l - p) + 2/i/_x(l - |) 

+/o(l-p)(l-/o) = ^(l-/o). (A5) 

Now, it is already known without assuming = 

[because the terms containing were canceled in 

Eq. Q] that in the ordered state fo = 2(1 —p)/p. Using 
this in the above equation and simplifying, we get 



= o. 



(A6) 



In the generalised version (Eq. ([3])) the time evolution 
equations for f\ and /_i will take the form 



dfi fofiq f2n w-. <7x , , u , pq, 

-£ = — - h (i -p)(i - 3) - - y) 



+(l-p)(l-|)}-/i/o(l-p), 



(A7) 



df-i fof-iq - 2 n „\ t -\ q^ £ £ r/-. PQ 



dt 



2 /^(i-pXi-i) -/,/_!{(! -if) 



+(l-p)(l-f)}-/-i/o(l-p). 



(A8) 



For steady state, the left hand sides will be zero. Then, 
as before, adding the resulting two equations and simpli- 
fying we get 

(1 - p)(l - |){(1 - / ) 2 - 2hf^} + 2/ 1 /_ 1 {(l - f ) 

+ (1 - |)} + /o(l - /o)(l -p) = ?f{l- MA9) 

Again, it was derived (without assuming polarization) 
that in the ordered state f = (p — l)(q — 2) /pq. Using 
this in the above equation and simplifying one would get 
/i/-i=0. 
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So we see that in the present model (and in its gen- 
eralised version), in the ordered state, ffj v _ 1 = for all 
x, y ^ as was used in the main text. Note that once 
one of the two extremes (—1 or +1) are completely elim- 
inated, it can not come back. So the non-coexistence is 
valid not only for 'space' but also for all times (after the 
steady state is reached). Note, however, that this con- 
dition (/1/-1 = 0) was assumed in Sec. IV and V and 



were numerically verified. 
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